Ground-state fidelity and quantum criticality in a two-leg ladder with cyclic four-spin exchange 
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We investigate a two-leg Heisenberg spin ladder with cyclic four-spin exchange by exploiting a newly- 
developed tensor network algorithm. The algorithm allows to efficiently compute the ground-state fidelity per 
lattice site, which enables us to establish the ground-state phase diagram for quantum lattice many-body sys- 
tems. The latter is based on the observation that, for an infinite-size system, any singularity on a ground-state 
fidelity surface characterizes a critical point, at which the system undergoes a phase transition. For the two-leg 
Heisenberg spin- 1/2 ladder with cyclic four-spin exchange, six different phases are identified: the ferromagnetic 
phase, the rung singlet phase, the staggered dimer phase, the scalar chirality phase, the dominant vector chirality 
region, and the dominant collinear spin region. Our findings are in a good agreement with the previous studies 
from the exact diagonalization and the density-matrix renormalization group. 
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Introduction. In the last decades, low-dimensional quan- 
tum spin systems, such as spin ladders, have been the sub- 
ject of extensive experimental and theoretical interest. Many 
fascinating features of the ladder systems have long been un- 
derstood theoretically from both analytical and numerical ap- 
proaches Uj,i2J- Among them is an intriguing property that 
the existence of an excitation gap depends on the number of 
legs: spin excitations are gapful for an even-leg ladder, and 
gapless for an odd-leg ladder There are also a number of low- 
dimensional cuprate compounds of transition metals, whose 
properties are described adequately by multi-leg spin ladders. 
These ladder systems often exhibit some attractive propeities 
such as quantum criticality. In addition, properties of several 
classes of materials such sls S rCu202 [3], La(,Caii,Cu2404i |0|, 
and {CsH\2N)2CuBr4, are well described by a two-leg Heisen- 
berg ladder with cyclic four-spin exchange. 

On the other hand, a novel approach to critical phenom- 
ena in quantum many-body physics has emerged, based on 
the fact that fidelity, a basic notion in quantum information 
science, is a measure of quantum state distinguishability. The 
approach allows us to characterize critical phenomena in a va- 
riety of quantum many-body lattice systems in any spatial di- 
mensions [5- 14]. As argued in Refs. ||6-11 1, the ground-state 
fidelity per lattice site is able to capture drastic changes of the 
ground-state wave functions around a critical point. This, in 
combination with the fact that many powerful numerical algo- 
rithms have been developed in the context of the tensor net- 
work (TN) representation, provides a powerful means to un- 
veil quantum criticality underlying quantum many-body sys- 
tems. In fact, a systematic scheme to study critical phenom- 
ena in quantum many-body lattice systems consists of three 
steps, as advocated in Ref. |15]: first, map out the ground- 
state phase diagram by computing the ground-state fidelity per 
lattice site; second, derive local order parameters (if any) from 
the reduced density matrices for a representative ground-state 
wave function in a given phase; third, characterize any phase 
without any long range order. 

In this paper, we consider a two-leg Heisenberg spin- 1/2 
ladder with cyclic four-spin exchange. The ladder system ex- 



hibits a very rich phase diagram, with six diff'erent phases 
identified: the feiTomagnetic phase, the rung singlet phase, 
the staggered dimer phase, the scalar chirality phase, the 
dominant vector chirality region, and the dominant collinear 
spin region, when a properly chosen control parameter is var- 
ied. Therefore, the ladder system provides a test bed for our 
scheme to study quantum criticality in spin ladders. This is 
achieved by exploiting a newly-developed TN algorithm lfl6ll . 
which allows to efficiently compute the ground-state fidelity 
per lattice site, and the reduced density matrix from a repre- 
sentative ground-state wave function. Our findings are in a 
good agreement with the previous studies from the exact di- 
agonalization ilTlllSl] and the density-matrix renormalization 
group (DMRgTusT 

Tensor network representation for spin ladders. The 
TN representation is a convenient way to represent ground- 
state wave functions in classical simulations of quantum 
many-body lattice systems, such as the matrix product state 
(MPS) Clil in one spatial dimension and the projected 
entangled-pair state (PEPS) ll24l - l27ll in two and higher spatial 
dimensions. Here, let us briefly recall the gradient TN algo- 
rithm to compute the ground-state wave functions for quan- 
tum many-body systems on an infinite-size two-leg spin lad- 
der 1 16], adapted to be suitable to a system with cyclic four- 
spin exchange. 

Assume that the Hamiltonian is translation-invariant: H - 
YAih^'\ with being the /-th plaquette Hamiltonian density 
along the leg direction. A TN representation for a quantum 
wave function consists of four-index tensors A 
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and Di attached to each site in the unit cell. Here, 5 is a 

tru ' 

physical index, s - \, d, with d being the dimension of the 
local Hilbert space, and £, r, u, and d denote bond indices, i, 
r, u, d - 1, D, with D being the bond dimension. Then, for 
a random initial state li/ro), the energy is a functional of the TN 
tensors: 

^ _ (MHl^o) 
(lAoltAo) 

which allows an efficient computation in the context of the TN 
representation. To update four-index tensors Ai ,, Bi,, C 
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and O^^.^^, we need to compute the energy gradient. 
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where X^^ e {Al^, B^^, D^J. Once the energy gradient 
is known, the four-index tensors A^^^, B^^^, C^^^^, and D^^.^^ may 
be updated as follows. 
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Here, 6 denotes the step size during each iteration, which 
is tuned to be decreasing in the implementation, when the 
ground-state wave function is approached. Here, we stress 
that four different tensors A^^^, B^^^, Cf.^^, and D^^^^ are up- 
dated simultaneously. Repeating this updating procedure until 
the ground-state energy converges, the system's ground-state 
wave function is generated in the TN representation. Actually, 
if the energy gradient |^ — > 0, a good approximation to the 
ground-state wave function is anticipated. 

The model. The two-leg Heisenberg spin- 1/2 ladder with 
cyclic four-spin exchange is described by the Hamiltonian: 



H - J± ^ S ij ■ S 2,i + J\\ ^^(S l,i ■ S ij+i + S 2,i ■ S 2,,+l) 



(4) 



Here, Sa,i (a = 1,2) denotes the spin- 1/2 Pauli operators at 
site / on the a-th leg, 7^ is the interchain coupling between two 
spins on each rung, is the intrachain coupling between two 
neighboring spins in each chain, and K is the cyclic four-spin 
exchange interaction coupling. The cyclic four-spin permu- 
tation operator exchanges the four spins around 
the !-th plaquette as Sij — > S ij+i — > 82,1+1 S2,i Sij, 
which can be decomposed in terms of the PauU spin operators 
involving bilinear and biquadratic terms: 

Pi,i+\ + Pij+\ -S 1,1 ■ S -I- S ij+i ■ S 2J+i + S 2,i+\ ■ S 2,i 
+ S 2,i ■ S ij + S ij ■ S 2,;+i + S ij+i ■ S 2,i 

+ 4(5i,, -5i,,+i)(52,,+i -52,/) (5) 

+ 4(S2,i-Sl,i)(Si,M-S2,M) 

- 4(5 1_/ ■ S 2,i+\)(S ■ S 2,i). 

The model has been inv estig ated b y the exact diagonaliza- 
tion 113128] and DMRG iHEllIiQ. Here, we focus on 
the computation of the ground-state fidelity per lattice site in 
terms of the newly-developed TN algorithm. For simplicity, 
we choose K - sin 6, and - J\\ = cos 6, with 6 e [-n, n]. 

The ground- state fidelity per lattice site. The ground-state 
wave function is generated from the newly-developed TN al- 
gorithm for a given choice of the coupling constants of the 
spin ladder, which allows to efliciently evaluate the ground- 
state fidelity per lattice site, a universal marker to detect quan- 
tum criticalities: a phase transition point is characterized by a 
pinch point on the fidelity surface. 



For the two-leg Heisenberg spin ladder with cyclic four- 
spin exchange, we choose as a control parameter For two 
different ground states, \^(9\)) and |i/'(02)), corresponding to 
two different values 9\ and 62 of the control parameter 6, the 
ground-state fidelity F(0i,02) - KiA(02)|iA(fi'i))l asymptotically 
scales as F(0i,02) ~ d{0\,02)^ , with the total number of 
the lattice sites. Here, d{0\,02) is the scaling parameter, in- 
troduced in Refs.|0-8[| for one-dimensional quantum lattice 
systems and in Refs.||9t] for two and higher-dimensional quan- 
tum lattice systems. In fact, d{0i,02) characterizes how fast 
the fidelity goes to zero when the thermodynamic limit is ap- 
proached. Physically, the scaling parameter d(0\,02) is the 
averaged fidelity per lattice site. 

In F{0i,02) 
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which is well defined in the thermodynamic limit even if 
F{0\,02) becomes trivially zero. It satisfies the properties in- 
herited from the fidelity F(0i,02): (i) normalization d(0, 0) = 
1; (ii) symmetry d{0\,02) - d{02,0\)\ and (iii) range < 

d{0x,02) < 1. 

In Fig[T] we plot a two-dimensional fidelity surface embed- 
ded in a three-dimensional Euclidean space for the two-leg 
Heisenberg spin-1/2 ladder with cyclic four-spin exchange. 
As shown in the upper panel, there are six pinch points on 
the fidelity surface, implying six phase transition points. In 
the lower panel, the contour plot of the ground-state fidelity 
per lattice site d{0\,02) on the 02)-plane is shown, with 
the truncation dimension up to 4. Therefore, we are able, by 
evaluating the ground-state fidelity per lattice site, to iden- 
tify six different phases: the ferromagnetic phase, the rung 
singlet phase, the staggered dimer phase, the scalar chirality 
phase, the dominant vector chirality region, and the dominant 
collinear spin region. Notice that, among six transition points, 
there are two first-order phase transitions at s; -0.40;7r 
and K 0.947r between ferromagnetic phase and its adjacent 
phases: the rung singlet phase and the dominant collinear spin 
region. The remaining four transition points are continuous. 
These results are in a good agreement with the earlier anal- 
yses 1 1/7, T3l based on the exact diagonalization and DMRG. 
Therefore, the TN algorithm yields reliable results for the two- 
leg Heisenberg spin ladder with cyclic four-spin exchange. In 
addition, the ground-state fidelity per lattice site, as a univer- 
sal marker to detect phase transitions, is able to capture drastic 
changes of ground-state wave functions around critical points 
for quantum spin ladders. 

Order parameters. Once the ground-state phase diagram 
is known, we are able to read out local order parameters (if 
any) from the reduced density matrices for a representative 
ground-state wave function in a given phase, as advocated in 
Ref. Ili. 

In the ferromagnetic phase, the non-zero-entry structure of 
the one-site reduced density matrix shows that the {Sa.i) are 
the same at all the lattice sites for the two-leg spin ladder 
Therefore, the local order parameter is 



Of = (lAol^Q-./liAo)- 



(7) 
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FIG. 1: (color online) Upper panel: A two-dimensional fidelity 
surface embedded in a three-dimensional Euclidean space. Lower 
panel: The countour plot of the ground-state fidelity per lattice site, 
d(Q[,6i), on the (6*1, 6'2)-plane, for the two-leg Heisenberg spin-1/2 
ladder with cyclic four-spin exchange. There are six pinch points 
on the fidelity surface. Therefore, six different phases are identified: 
the ferromagnetic phase, the rung singlet phase, the staggered dimer 
phase, the scalar chirality phase, the dominant vector chirality re- 
gion, and the dominant collinear spin region. The ground-state phase 
diagram is as follows: the ferromagnetic phase for -1.06;n'(0.94;n') < 
9 < -OAQn, the rung singlet phase for -OAOn < 9 < 0.06;r, the 
staggered dimer phase for O.Obn < 9 < 0.15;r, the scalar chirality 
phase for O.lSn < 9 < 0.38;r, the dominant vector chirality region 
for 0.38;r < 9 < 0.80;r, and the dominant collinear spin region for 
0.80;r <9< 0.94n. 



As seen in Fig.|2]l, \0f\ = 0.50. That is, spins are fully polar- 
ized in this phase. In fact, the spin correlations {Saj ■ Saj+\} 
between the nearest-neighbor spins on the legs and the spin 
correlations (S ij ■ 82,1) between spins on the rungs are 0.25. 
The ferromagnetic state minimizes the energy on each plaque- 
tte separately for -1.067r(0.947r) <d< -0.407r. 

In the rung singlet phase, the ground-state wave function 
may be approximated by the product of local rung singlets. 
The ground-state lacks long-range order in the conventional 
sense, thus there is no local order parameter; instead, an exotic 
order occurs. In fact, the rung singlet phase and the Haldane 
phase are essentially the same, namely, both are characterized 
by the so-called string order ||32| . The rung singlet phase lies 
in -0.407r <0< OMn. 

In the staggered dimer phase, the non-zero-entry structure 



of the two-site reduced density matrix exhibits a pattern, with 
the local order parameter as follows, 

OsD =r-(i/'o|5 ■ S ij - S ij ■ S 

Z (8) 

+ S 2.i ■ S 2,/+l - S 2J-1 ■ S 2,/|l/'o)- 

Here, {i//o\Sa,i ■ Saj+ilif/o) = (lAol^ 1,1 • S i.iliAo) if a + i is even, 
{tl/o\S a.i-S aj+i\^o) = (i/'o|52,r52,2|i/'o)if<3'+'is oddfor two de- 
generate symmetry-breaking ground states in this phase. The 
ladder is in the staggered dimer phase for 0.067r < < 0.157r 
(cf. FiglDi. 

In the scalar chirality phase, we need to study the non-zero- 
entry structure of the three-site reduced density matrix. This 
yields the local order parameter 

Osc = ('Ao|5i,,-(52,,x5i,,+i)|(Ao). (9) 

It breaks the spatial symmetries and the time reversal symme- 
try, but not the internal SU(2) symmetry. The scalar chirality 
phase lies in 0.157r < 6 < 0.387r, as seen from Fig|2] 

In the dominant vector chirality region, the non-zero-entry 
structure of the two-site reduced density matrix yields the lo- 
cal order parameter 

Ovc-leg = (tAol^a,/ X S a,i+l\l^o), (10) 
Ovc-rung = {lAol^a./ X S a+lMo)- (H) 

It breaks the spatial symmetries and the time reversal sym- 
metry. In this phase, the spin correlations are strong between 
bonds on rungs and legs, but the spin correlations are very 
weak between diagonal bonds. The order parameter is plotted 
in Figlll which is non-zero between 0.38?: < 6* < O.SOtt. 

In the dominant collinear spin region, spins on the same leg 
exhibit ferromagnetic correlations, while spins on the same 
rung exhibit antiferromagnetic correlations. The non-zero- 
entry structure of the one-site reduced density matrix yields 
the local order parameter 

Ocs =^{MSi,i-S2MQ)- (12) 

The dominant collinear spin region lies in 0.80;7r ^ ^ 0.94.7r. 

Therefore, we are able to "derive", by investigating the non- 
zero-entry structure of the reduced density matrices for repre- 
sentative ground-state wave functions from different phases, 
the local order parameters for the ferromagnetic phase, the 
staggered dimer phase, the scalar chirality phase, the domi- 
nant vector chirality region and the dominant collinear spin 
region, with the order parameter Of, Osd, Osc, Ovc and Ocs 
explicitly shown in Fig|2l In addition, no local order param- 
eter is detected in the rung singlet phase, indicating that long 
range order is lacking in this phase. The ground-state phase 
diagram established from the local order parameters coincides 
with that from the ground-state fidelity per lattice site. That is, 
the ladder system undergoes four continuous phase transitions 
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FIG. 2: (color online) The local order parameters Of, Osd, Osc, 
Ovc, and Ocs in the ferromagnetic phase , the staggered dimer phase, 
the scalar chirality phase, the dominant vector chirality region, and 
the dominant coUinear spin region versus 9, respectively. 



at * 0.06;r, ^ O.IStt, 6 ^ 0.38;r and 6 ~ O.SOtt, and two 
first-order phase transitions at ^ -0.40?: and 6 » 0.947r. 

Conclusions. We have exploited a newly-developed effi- 
cient TN algorithm to compute ground-state wave functions 
for the two-leg Heisenberg ladder with cyclic four-spin ex- 
change. We have mapped out, by computing the ground-state 
fidelity per lattice site, the phase diagram of the two-leg spin 
ladder with cyclic four-spin exchange. Six different phases are 
identified: the ferromagnetic phase, the rung singlet phase, the 
staggered dimer phase, the scalar chirality phase, the domi- 
nant vector chirality phase, and the dominant coUinear spin 
phase. In addition, the corresponding local order parameters 
are "derived" from the reduced density matrices, which are 
computed efficiently in the context of the TN representation. 
Our findings are in a good agreement with the previous stud- 
ies from the exact diagonahzation 113, Ell and DMRG jisll . 
Therefore, the TN algorithm yields reliable results, and the 
ground-state fidelity per lattice site, as a universal marker to 
detect phase transitions, is able to capture drastic changes of 
ground-state wave functions around critical points, even for 
the two-leg Heisenberg spin ladder with cyclic four-spin ex- 
change. 
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